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Abstract. Deep images of the galactic globular cluster 
M 4 taken at various locations with the NIC 3 and the 
WFPC 2 cameras on HST were used to derive detailed lo- 
cal optical and near IR luminosity functions. White dwarfs 
have been detected for the first time on a color sequence at 
constant luminosity in the F110W band. Transformation 
of the observed luminosity functions into mass functions 
via the most up to date theoretical mass luminosity rela- 
tions currently available results in best fit local mass func- 
tions, in logarithmic mass units, that consist of a power- 
law dN oc m~ x d\ogm, with single exponent x = —0.8 for 
the inner regions and a two-segment power-law which rises 
with x = 1 down to Msi4 ~ 8.5 and then drops all the way 
to the detection limit with x — —0.3 for the outer regions. 
This behavior cannot be reconciled with the expectations 
of a multi-mass King-Michie dynamical model using as 
input the canonical structure parameters for this cluster 
(core radius r c = 50" and concentration c = 1.6; Harris 
1996). Thus, either the model does not accurately reflect 
the structure of the cluster due to some effect not prop- 
erly accounted for in it or the canonical cluster structural 
parameters have to be significantly modified. Reasonable 
fits to all the present observations can be obtained with 
various global mass functions provided the cluster's struc- 
tural parameters such as concentration c and core radius 
r c are in the range c G [1.4,1.9] and r c <E [50", 90"]. 
The best compromise, in this case, consists in a model 
with a two-segment power-law mass function with expo- 
nents x = 0.2-1.0 in the mass interval 0.8 < M Q < 0.25, 
x = —0.4 for m < 0.25M Q and structural parameters that 
require the least modification from the currently estab- 
lished values. This last result differs only minimally from 
that obtained for other globular clusters studied so far 
with HST which seem to have global mass functions in- 



Send offprint requests to: F. Paresce 

* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 
which is operated by AURA for NASA under contract NAS5- 
26555 

** On leave from Osservatorio Astronomico di Roma, Via 
dell'Osservatorio 2, 1-00040 Monte Porzio Catone (Roma), 
Italy 



creasing up to a peak at ~ 0.25M Q and then flattening 
out and possibly dropping to the H-burning limit. 

Key words: stars: Hertzsprung-Russel (HR) and C-M 
diagrams - stars: luminosity function, mass function - 
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1. Introduction 

The galactic globular cluster M4 (NGC6121) is well 
suited to detailed measurements of its faint stellar pop- 
ulation due to its proximity (~ 2kpc) and looseness 
(c = 1.59; Harris 1996) that allows one, in principle, to 
probe its main sequence (MS) down to close to the H- 
burning limit with relative ease even from the ground 
(Kanatas et al. 1995). This capability is crucial in deter- 
mining the shape of the stellar initial mass function (IMF) 
at the low mass end where observations have been noto- 
riously sparse and unreliable (Chabrier & Mera 1997, De 
Marchi & Paresce 1997). In particular, it is this part of 
the mass function (MF) close to the peak of the luminos- 
ity function that carries the most critical information con- 
cerning the possible log-normal form of the IMF (Adams 
& Fatuzzo 1996, Scalo 1998). Moreover, it is also here 
that one expects the effects of mass segregation and tidal 
stripping to be most severe and, therefore, most visible. 
If measurable, this signature could tell us much about the 
extent and nature of this phenomenon which is only im- 
perfectly understood at the moment (Gnedin & Ostriker 
1997, Vesperini & Heggie 1997). 

In practice, access to this region of the MS is par- 
ticularly difficult due to a number of factors affecting the 
precision with which a luminosity function (LF) can be de- 
termined and a MF derived from the observed LF. Incom- 
pleteness, field star contamination, small number statis- 
tics, imperfect calibrations and, especially, very uncertain 
mass-luminosity relations (M-L) for low-mass stars have 
so far all conspired to severely limit the usefulness of M 4 
for this purpose. Kanatas et al. (1995) have come tantaliz- 
ingly close to the end of the MS of M 4 but the reliability 
of their results has been difficult to assess owing especially 
to the crudeness of the M-L relationship and to the inabil- 
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ity to correct observationally for field star contamination. 
The most interesting aspect of their work has been the sug- 
gestion that this cluster is particularly poor in low mass 
stars which might imply that it has been significantly dis- 
rupted by internal and/or external forces. IR photometry 
from the ground carried out by Davidge & Simons (1994) 
has also found a flat MF down to O.2M . 

In view of the recent availability of precise inter- 
nally consistent M-L relationships for these metallicities 
(Baraffc et al. 1997) and the overwhelming importance to 
have at least one other cluster to compare with NGC 6397, 
the only one for which, up to now, we had a reliable MF 
extending all the way down to the H-burning limit and 
to establish clearly the signature of tidal stripping on any 
cluster, we decided to observe M4 with the NICMOS in- 
strument on board HST and to complement this exten- 
sive data set with archival WFPC 2 observations in order 
to achieve a double goal: extend the LF reliably down to 
the H-burning limit and cover as wide an area in the clus- 
ter as possible to correctly account for the effects of mass 
segregation. In this paper, we present the results of this 
investigation. Details of the observations, their analysis 
and the color-magnitude diagrams are described in Sect. 2 
in which we compare the NICMOS observations with the 
WFPC 2 results, showing how the infrared photometry has 
reached MS stars close to the H-burning limit, and we con- 
firm the flattening of the MF for masses m < 0.25 M© . 
The LFs and MFs are examined in Sect. 3 whereas Sect. 4 
shows how the various LFs have been used to constrain a 
multi-mass King-Michie dynamical model of M 4. A sum- 
mary of the results follows in Sect. 5. 

2. Observations and analysis 

The regions of the cluster analyzed in this paper are shown 
in Fig. 1. The two HST-NIC3 images (dashed squares) 
were taken at ~ I'. 9 and ~ 2f9NE of the center. The 
two HST- WFPC 2 archive fields used here were located 
~ 2' and ~ 6' E of the cluster center ("L" shaped boxes). 
In the following, we describe these observations in detail. 

2.1. The NICMOS data 

The observations of M 4 described here were obtained on 
1998 January 31 and February 1 (UT) with the NICMOS 
instrument (MacKenty et al. 1997) on board the HST, 
during the first so called "NIC 3 Campaign," when the 
HST secondary mirror was moved to place the NIC 3 cam- 
era into focus. All the images were taken with this camera, 
which has a plate scale of O'.'2/pix, yielding a field of view 
of 51'/2 x 5172. 

Four exposures were taken at RA = 16 : 23 : 42.6 
and DEC = -26 : 30 : 40, or - L'9NE of the cluster 
center (field FN 1). Four more exposures were taken, with 
the same instrumental setup, at RA — 16 : 23 : 47.6, 
DEC = -26 : 30 : 30, or ~ 2f9NE of the cluster center 
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Fig. 1. Location of all the fields used in this paper in M4 
with respect to the cluster center marked with a cross. 
The dashed radial circles start at 1' with steps of 1'. The 
dashed squares are the 51" x 51" wide NIC3 fields and 
the two fields bounded by the thin lines were observed 
through the WFPC 2 

Table 1. NIC 3 observations of the cluster M4 



field 


filter 


texp 


Sampl.Seq. 


N. Sampl. 


N. exp 


FN 1 


F110W 


1919.87 


STEP 128 


24 


x2 


FN 1 


F160W 


703.94 


STEP 28 


20 




FN 1 


F160W 


767.94 


STEP 28 


24 




FN 2 


F110W 


1919.87 


STEP 128 


24 


x2 


FN 2 


F160W 


767.94 


STEP 28 


24 


x2 


SKY 


F110W 


1919.87 


STEP 128 


24 


x2 


SKY 


F160W 


703.94 


STEP 28 


20 




SKY 


F160W 


767.94 


STEP 28 


24 





(field FN 2). Finally, two additional images were taken at 
RA = 16 : 24 : 04.7, DEC = -26 : 36 : 26 (SKY field), i.e. 
2° away from the cluster center to account for the contam- 
ination due to field stars. All these images were obtained 
through the F110W and F160W filters using the MULTI- 
ACCUM mode. Table 1 summarizes the main characteris- 
tics of the images of these three fields. For each observed 
field as indicated in column 1, we report the adopted fil- 
ter, the exposure time in seconds, the MULTIACCUM 
predefined sample sequence, the number of exposures for 
each sample sequence and the number of MULTIACCUM 
sequences for each field. 

All frames have been processed with the standard 
NICMOS-HST pipeline procedure (CALNICA) to per- 
form bias subtraction, linearity correction, dark count cor- 
rection, flat fielding and cosmic ray identification and re- 
jection. The frames in the same field and taken through 
the same filter have been registered with respect to each 
other and averaged to improve the statistics. The final 
images correspond to a total exposure time of 64 min for 
the F110W filter in all the observed fields, 24 min for the 
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Fig. 2. Negative image of the field FN 1 in the globular 
cluster M4, located ~ If 9 away from the nominal cluster 
center, as observed with the F110W filter of the NIC 3 
camera. The field covers ~ 51" on a side. North is in- 
dicated by the arrow in the upper-left corner with East 
indicated to the left of North. The exposure time of the 
combined image is 3,840 s 

F160W filter in FN 1 and SKY and 27 min for the F160W 
filter in FN 2. 

The FN1, FN 2 and SKY combined frames are shown 
in Figs. 2, 3, and 4, respectively. The image quality is excel- 
lent and the background is almost flat with the exception 
of the regions surrounding bright stars. The average per 
pixel count rate is 0.030 count s _1 and 0.035 count s^ 1 , re- 
spectively in F110W and F160W, with standard deviation 
of 0.015 count s _1 and 0.011 count s^ 1 . 

Photometry was carried out using the standard IRAF 
DAOPHOT package. The automated star detection rou- 
tine daophot.daofind was applied to the frames, by setting 
the detection threshold at 5 a above the local background 
level. We have then examined by eye each individual ob- 
ject detected by daofind and discarded a number of fea- 
tures (PSF tendrils, noise spikes, etc.) that daofind had 
misinterpreted as stars, as well as a few extended objects 
whose full width at half maximum (FWHM) exceeded by a 
factor of two or more that typical of stellar objects in our 
frames (FWHM ~ 1.3 pixel). As regards the saturation 
of bright objects, the CALNICA pipeline verifies whether 
and when the signal in each pixel, in the multiple read- 
outs of each MULTIACCUM sequence, reaches the satura- 
tion or nonlinearity limit prior to the end of the exposure. 
All the samples which occur after the saturation limit are 
identified as bad and are not used in the final calculation 
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Fig. 3. Same as Fig. 2 but for Field FN 2, located 2'. 9 away 
from the center of the cluster 
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Fig. 4. Same as Fig. 2 but for Field SKY 



of the mean count rate. This procedure should, in princi- 
ple, allow one to measure fluxes in a wide dynamic range 
(~ 10 mag in our case). 

A total of 466, 332 and 175 objects respectively in the 
FN 1, FN 2 and SKY fields were identified and their fluxes 
measured with the PSF-fitting daophot. allstar photometry 
package, to derive the color-magnitude diagrams (CMD) 
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for the three fields. The internal photometric accuracy has 
been estimated both with the IRAF phot task and through 
standard artificial-star experiments, giving similar results: 
better than ±0.02 mag for bright stars, degrading to about 
±0.3 mag toward the detection limit at mno ~ 27. 

Instrumental magnitudes were calibrated and con- 
verted into the HST magnitude system (STMAG) using 
the relation: 

m ST = -2.5 log (^pj -21.1 (1) 

where c is the count rate measured for each star, U the 
inverse sensitivity of the instrumental setup (camera + 
filters), and s the encircled energy (i.e. the fraction of 
the total flux sampled by our photometry). In order to 
estimate the encircled energy e, we have used five iso- 
lated, relatively bright stars, and have compared the flux 
measured by our photometry to their total flux, defined 
as that falling within a 2" radius. We have found in this 
way that, for our choice of aperture radius (2.5 pixel) and 
background annulus (2.5-10 pixel) the encircled energy 
amounts to ~ 90% for both the F110W and F160W fil- 
ters. Because the calibration of the NIC 3 camera is still 
preliminary, and the values of the inverse sensitivity U 
are based on a limited set of spectrophotometric stan- 
dards, we cannot convert very accurately our magnitudes 
from the STMAG system into any ground-based system. 
We, therefore, expect our absolute photometry to be ac- 
curate to within 10-15% (Colina & Rieke 1997). It is still 
possible, however, to translate our measurements into the 
VEGAMAG photometric system of the HST, defined as 
one in which the magnitude of Vega would be in all bands 
and which is more similar to the classical Johnson-Cousin 
ground-based system. This can be done by simply sub- 
tracting the zero-point constants of 2.3 mag and 3.7 mag 
from the values of mno and mieo, respectively, thus ob- 
taining the J and H magnitudes that we use hereafter. 

The CMDs of all the objects identified in the three 
fields observed (FN 1, FN 2 and SKY) are shown in Figs. 5, 
6, and 7. The MS is well defined in both cluster fields 
down to J ~ 20, but below this level the stellar density 
decreases abruptly. In the SKY field, located 2° away from 
the cluster center, the density of stars begins to increase 
for J > 20, implying that the stellar population of the 
cluster becomes negligible at this point. Using the M-L 
relation of Baraffe et al. (1997), this point corresponds to 
~ 0.2 M Q for a distance modulus of (m— M)j = 11.86. An- 
other noticeable feature of the MS is the change of slope at 
J ~ 18 corresponding to M ~ 0.5M Q . This is because of 
the increasing infrared absorption of H2 molecules in the 
stellar atmosphere which counteracts the reddening due to 
both the decreasing temperature and increasing metallic 
molecular absorption in the optical band (Baraffe et al. 
1997). Below this limit, the MS becomes almost vertical. 

The broadening of the MS for J < 16.5 is greater than 
expected on the basis of the photometric errors alone. 



This broadening has been already noticed by Richer et 
al. (1995). Although it cannot be excluded that this fea- 
ture is partly due to differential reddening (Alcaino et al. 
1997) or to the presence of a binary sequence (Kroupa & 
Tout 1992; Ferraro et al. 1997), it seems more likely that 
this effect has to be ascribed to the non uniformity in the 
pixel response function of the camera (intra-pixel sensi- 
tivity) which has been shown to affect the photometry at 
the 10% level in NIC 3 (Storrs 1998). On the other hand, 
the CMD of M4 published by Kanatas et al. (1995) does 
not show this peculiar MS broadening. We have estimated 
that if the MS spread were to depend on the data process- 
ing alone, the actual error should be almost four times 
larger than the photometric uncertainty (~ 0.04 mag) in 
the range J < 16.5. As a consequence, it would be difficult 
to evaluate the age of the cluster from these data, but the 
LF of MS stars is little affected by this systematic effect. 

We have superposed on the data shown in Figs. 5 and 
6 the expected lOGyr, [M/H]= —1.0 theoretical isochrone 
(Chabrier & Baraffe 1997; Allard & Hauschildt 1997), as 
seen through the NICMOS filters F110W and F160W. For 
this comparison we have adopted for M4 a color excess 
E(B-V)= 0.40 and a distance modulus of (m — M) Q = 
11.51 (Djorgovski & Meylan 1993). The ensuing values 
of the absorption in the F110W and F160W bands, es- 
timated with the relation of Cardelli, Clayton, & Mathis 
(1989), are respectively A 110 = 0.33 and A 160 = 0.18. Al- 
though these values do not bring the theoretical models 
into perfect agreement with the observed CMD (an addi- 
tional ~ 0.1 mag shift to the red seems to be required), 
we regard this mismatch as insignificant in light of the 
uncertainties on both the cluster's parameters and the- 
ory. Moreover, there are indications of higher and variable 
values for Ry (Vrba, Coyne, & Rapia 1993) towards the 
direction of M 4 which is behind the Sco-Oph dust com- 
plex. Thus, because the calibration of our filters is still 
preliminary, the match between observation and theory 
as shown in both the FN 1 and FN 2 fields (Figs. 5 and 6) 
is remarkable. 

The completeness of the sample has been evaluated by 
running artificial star tests in both bands. For each 0.5 
magnitude bin we have carried out 10 trials by adding a 
number of artificial stars not exceeding 10 % of the total 
number of objects in that bin and using a PSF derived di- 
rectly from the co-added frames. These trials were followed 
recursively by daofind and allstar runs with the same pa- 
rameters used in the reduction of the scientific images to 
assess the fraction of objects recovered by the procedure 
and the associated photometric error. Moreover, for each 
artificially added star, we compared the input magnitude 
with that at which the star was recovered, and found no 
systematic difference. In addition, the standard deviation 
of the magnitude values at which stars are recovered is in 
excellent agreement with the 1 a uncertainty of the pho- 
tometric reduction, as estimated by the IRAF phot task, 
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Fig. 5. Color-magnitude diagram of the field FN 1 located 
l! 9 off the center of the globular cluster M4 (NGC 6121). 
A total of 466 stars are detected and measured in two 
bandpasses (F110W and F160W) with the NIC 3 camera 
in this field. The solid line represents the expected theo- 
retical MS from Baraffe et al. (1997). The magnitudes and 
colors in the HST-VEGAMAG system are shown respec- 
tively on the right and on the upper axis of the diagram. 
Photometric error bars are indicated on the left. The stel- 
lar masses determined by the M-L relationship of Baraffe 
et al. (1997) are indicated on the right side of the diagram 



at all magnitude levels and is indicated on the left-hand 
side of the CMDs shown in Figs. 5, 6, 7. 

The photometric completeness of the three observed 
fields is displayed in Fig. 8 as a function of the rano mag- 
nitude, whereas Table 2 lists, for each magnitude bin, the 
completeness fraction, the corrected LF, and the corre- 
sponding rms errors coming from the poissonian statistics 
of the counting process (all values of the LFs as well as 
their errors have been rounded off to the nearest integer). 

Because we planned the NICMOS exposures in such 
a way that the F110W and F160W band images reach 
the same limiting magnitude with the same SNR, and be- 
cause the colors of all the objects seen in the CMD (both 
cluster and field stars) span a very narrow range, all the 
objects identified in the F110W frames are also found in 
the F160W images and, therefore, the photometric com- 
pleteness is similar in both bands. 

To determine the stellar LF of M 4 at the locations cor- 
responding to our fields FN 1 and FN 2 we have used the 
color and magnitude information derived from the CMDs 
in Figs. 5 and 6. Because M4 is projected towards the 
Galactic bulge at a low Galactic latitude (b ~ 16), the 



Fig. 6. Same as Fig. 3 but for field FN 2, located 2'. 9 away 
from the center of the cluster. A total of 332 objects are 
measured in this field 
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Fig. 7. CMD of the 175 objects in the field SKY at - 2° 
from the center of the cluster 



contamination due to field stars in FN 1 and FN 2 is ex- 
pected to be significant. While at optical wavelengths it 
is often possible to disentangle the stars belonging to the 
cluster from those that do not simply on the basis of their 
color (see next section for a description of the "2.5 a clip- 
ping" criterion), the colors of low-mass stars in the near-IR 
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83 


55 
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91 
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98 


11 


3 


22.25 


82 


29 


6 


90 


16 


4 


97 


8 


3 


22.75 


78 


27 


6 


88 


16 


4 


97 


13 


4 


23.25 


74 


22 


6 


85 


15 


4 


96 


17 


4 


23.75 


70 


17 


5 


81 


21 


5 


94 


13 


4 


24.25 


66 


26 


7 


75 


23 


6 


93 
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58 
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Fig. 8. Photometric completeness for the three observed 
fields. The completeness was calculated as the ratio of the 
numbers of the recovered stars N r and the added artificial 
stars N a 

make this a very difficult undertaking. And indeed, both 
cluster and field stars share the same narrow color region 
in our F110W, F160W CMDs (see Figs. 5, and 6), with 
the latter scattered in luminosity because of their vary- 
ing mass as well as distance. We had, however, expected 
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Fig. 9. The luminosity functions corrected for photomet- 
ric incompleteness in the three observed fields, FN 1 (solid 
line), FN 2 (dashed line) and SKY (dash-dotted line) 



this behavior (see Baraffe et al. 1997), and had therefore 
obtained a comparison field (SKY) at the same projected 
distance from the Galactic bulge as the FN 1 and FN 2 
fields. Thus, the use of a control field is mandatory to re- 
liably measure the stellar LF of this cluster, due to its 
unfavorable location. 

We have measured the number of stars in each 0.5 
magnitude bin as a function of the mno magnitude along 
the MS in all three CMDs, discarding stars bluer than 
J — H ~ 0.4 as they are clearly not MS objects. We 
have adopted exactly the same color selection also for the 
SKY field, so as to ensure the consistency of the compar- 
ison with the FN 1 and FN 2 fields data. Fig. 9 shows the 
LFs in the three observed fields FNl (continuous line), 
FN 2 (dashed line) and SKY field (dot-dashed line) after 
correction for photometric incompleteness as reported in 
Fig. 8 and Table 2. The artificial star tests indicate that 
at J ~ 23 on average more than 50%, 60%, and 80% 
of the added stars are recovered respectively in the FN 1 , 
FN 2, and SKY, thus implying statistically reliable LFs 
towards the H-burning limit at ~ O.IMq, corresponding 
to J = 22.5 according to the M-L relations of Baraffe 
et al. (1997). In spite of the different crowding level and 
photometric completeness, the LFs measured in the ad- 
jacent fields FN 1 and FN 2 match each other very well, 
both showing peaks at J ~ 16.5 and J ~ 19.5. It is inter- 
esting to note that, after correction for incompleteness, all 
three LFs tend to a plateau for J > 21, clearly marking 
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Fig. 10. IR luminosity functions of M4, obtained from our 
photometry and corrected for both incompleteness and 
contamination. The solid histogram represents the inner 
field (FN 1) LF whereas the dashed histogram reflects the 
outer field (FN 2). 

in this way the bottom of the MS, from where field star 
contamination increases. 

The final step in obtaining the LF of M 4 at these lo- 
cations in the cluster is to subtract, after correction for 
photometric incompleteness, the SKY LF from the LFs 
measured in FN 1 and FN 2. The LF derived in this way 
is shown in Fig. 10 where the LFs in fields FN 1 (solid 
line) and FN 2 (dashed line) are corrected for incomplete- 
ness and contamination. The more external LF (FN 2) has 
been shifted vertically to reach the same maximum as the 
other. Within the observational errors, the two LFs are 
indistinguishable from one another. 

2.2. The WFPC2 data 

HST WFPC 2 observations of M 4 are available in the 
HST archive, from which we have retrieved two sets of 
F555W, F814W band images obtained in 1995 at ~ lf3E 
(field FW 1) and - 5' SE (field FW 2) of the cluster center 
(see Fig. 1). Since field FW 1 contains our NICMOS FN 1 
field, we are able to compare the LFs measured at the 
same location in different bands. In addition, field FW2 
provides information on the variation of the LF as a func- 
tion of the radial distance. 

The raw images were processed using the standard 
HST pipeline calibration. Images of the same field taken 
through the same filter were combined to remove cosmic 
ray hits and to improve the SNR of the data. The total 



exposure time of the combined images of field FW 1 corre- 
sponds to 15, 000 s in the F555W band and 5, 500 s in the 
F814W band. Exposure times for FW2 arc 31, 500 s and 
7, 200 s respectively in F555W and F814W. 

We have measured the fluxes of all the unresolved ob- 
jects detected in these images by using almost the same 
photometric reduction technique employed for the NIC- 
MOS data (see Sect. 2.1), excluding the PC chips. In par- 
ticular, we have run the automated star-detection routine 
daofind on the combined frames to detect objects rising 
at least 5 a above the local average background. Each of 
the detected objects has been visually examined and the 
spurious identifications, due essentially to the diffraction 
spikes of a number of saturated stars, were removed from 
the output list. In this way a total of 1439 and 1078 unsat- 
urated objects were detected, respectively in fields FW 1 
and FW2. 

The photometric routines phot and allstar were used 
to measure the fluxes of the objects detected in such a way 
on the single frames. The PSF was modeled by averaging 
moderately bright stars for each frame. The instrumental 
magnitudes from the different frames were weighted by 
the inverse square of their uncertainties, averaged for each 
filter and converted to the instrumental magnitudes of the 
combined images. 

We calculated the aperture correction for each WF 
chip in the two fields, using a diameter size of 1" and 
transformed the final instrumental magnitudes into the 
VEGAMAG photometric system adopting the updated 
zero points listed in the November 1997 edition of the 
HST Data Handbook. 

Fig. 11 displays the msi4, m^ — msn CMDs for all ob- 
jects of the inner field FW 1 and of the external field FW 2. 
The completeness corrections have been estimated by run- 
ning artificial star tests in both bands and in the combined 
images for each WF chip. The procedure is essentially the 
same as described in Sect. 2.1. Our simulations indicate 
that in the inner field at msu — 21 on average about 90% 
of the artificial stars are recovered. In the outer field, the 
same completeness level is reached at m 8 i 4 ~ 22 due to 
the lower crowding. 

Fig. 12 shows the NICMOS and the WFPC 2 CMDs of 
the 193 stars common to both inner fields. The center of 
the overlapping region is located ~ \[ 9 NE of the cluster 
center. The crosses represent objects already identified as 
white dwarfs (WD) by Richer et al. (1997) on the basis 
of the agreement between observational and theoretical 
loci of a 0.5M Q Hydrogen-rich WD cooling sequence, the 
asterisks indicate objects which, lying outside the narrow 
cluster MS, probably belong to the foreground whereas 
the circles mark the cluster members. 

As a first result we have the observational confirmation 
of the contamination of the cluster MS in the infrared due 
to field stars: it is impossible to follow the MS in the IR 
bands all the way to its very bottom without a statistical 
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decontamination through the use of an external field or 
excluding the foreground stars using proper motion data. 

Secondly, the cross-identification of the stars in both 
the WFPC 2 and NICMOS photometry tables enables us 
to identify in a globular cluster for the first time an in- 
frared WD sequence. In the F110W band the WDs are 
distributed in a color sequence instead of a luminosity se- 
quence as in the optical CMD. This is an important obser- 
vational constraint for the modeling of the atmosphere of 
cool WDs. Bergeron, Saumon, & Wesemael (1995) have 
shown that the cool WDs with mixed H/He chemical 
composition can be easily recognized from their predicted 
strong infrared flux deficiency. 

Adopting the 2.5 a clipping criterion (De Marchi & 
Paresce 1995), from the CMDs of Fig. 11 we have mea- 
sured the LF of MS stars by counting the objects in 
each 0.5 mag bin and within ±2.5 times the 771555 — m 8 i4 
color standard deviation around the MS ridge line. The 
WFPC 2 LFs of both observed fields, corrected for in- 
completeness and field star contamination, are shown in 
Fig. 13. The exposure times of the WFPC2 frames were 
calculate to reach the very end of the MS and, as a conse- 
quence, the saturation effect does not allow us to evaluate 
reliable LFs at m 8 i4 < 19 while the foreground contami- 
nation cuts them at m^u > 22.5. 

3. The luminosity and mass functions 

Fig. 14 shows the F814W LFs of the two WFPC 2 fields 
located at ~ 2' (triangles) and at ~ 6' (boxes) from the 
cluster center, together with the NICMOS LF (diamonds) 
converted into the F814W band using the Baraffe et al. 
(1997) M-L relationships appropriate for the metallicity 
of M4 ([M/H] = -1, i.e. [Fe/H] ~ -1.3; Harris 1996) as 
a function of the absolute M 8 i4 magnitude. Through the 
theoretical M-L relations, we can associate a mass to each 
of the magnitude bins that appear in the NICMOS LF. 
Then, by also knowing the M-L relation for the F814W 
band, we can easily translate the center of the magnitude 
bins from the F110W magnitude scale to the F814W band 
and adjust the width of the bins themselves to reflect the 
different slope of the M-L relations in the magnitude range 
covered by each bin. 

Within the error bars, in the overlapping range, the 
inner NICMOS and inner WFPC 2 LFs show the same 
trend: the overall LF is relatively flat or slightly decreasing 
in the range 4 < M S u < 8 and decreases steeply beyond 
that point. The agreement between the two LFs measured 
in the visible (F814W) and in the near IR (F110W), in 
spite of the small size of the observed sample, confirms the 
robustness of our photometric reduction and the reliability 
of the M-L relations used to translate the IR data to the 
visible (see also De Marchi 1999 on this issue). 

On the other hand, the WFPC 2 LF measured 6' away 
from the center (~ 1.6rh m ; Harris 1996) shows a different 
shape, in that it increases with decreasing luminosity up 



2.5 




1.0 1 1 1 1 

IB 19 20 21 22 23 
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Fig. 14. Logarithm of the WFPC2 2' (triangles) and 6' 
(boxes) LFs in the F814W band together with the NIC- 
MOS LF (diamonds) translated to F814W adopting M-L 
relations of Baraffe et al. (1997). Solid lines represent the 
fit obtained using a power-law with exponent x = —0.8 
in the inner field. As regards the outer field, the LF is re- 
produced using a two-segment power-law, with exponents 
x = 1 down to M 8 i4 ~ 8.5 and x = —0.3 down to the 
detection limit 

to a peak at M 8 i4 ~ 9 and from there it drops all the 
way to the MS detection limit at M 8 i 4 ~ 10.5. This be- 
havior is typical of all the GCs observed so far with the 
HST at or near the half-mass radius (Paresce, De Marchi 
& Romaniello 1995; De Marchi & Paresce 1995a, 1995b, 
1996, 1997; Cool, Piotto, & King 1996; King et al. 1997; 
Pulone, et al. 1998, De Marchi 1999) and is related to an 
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extremum in the derivative of the mass-absolute magni- 
tude relation of the low mass stars as discussed in Kroupa 
& Tout (1997). 

Thus, the obvious difference between the inner and 
outer LFs should be compatible with the dynamical mod- 
ifications of a GC as a result of mass segregation due to 
energy equipartition through two-body relaxation (Spitzcr 
1987). Indeed, it is expected that more massive stars give 
up part of their kinetic energy to lighter stars during close 
encounters and sink towards the cluster center. This ef- 
fect would explain qualitatively why massive stars are rel- 
atively more numerous at 2' from the center than they are 
near the half-mass radius. We address this issue in more 
detail in Sect. 4. 

Under the assumption that the MF of the stellar pop- 
ulation at 2' from the cluster center is represented by an 
exponential distribution of the form: 

dN oc m~ x dlogm (2) 

and using the M-L relation of Baraffe et al., we can re- 
produce rather accurately the inner LF in Fig. 44 with 
an exponent x = —0.8 ± 0.4 (solid line). The outer LF, 
however, cannot be fit by any single exponent power-law 
distribution, regardless of the adopted slope, as no such 
function can reproduce the peak observed at M^u — 9. 
A reasonable fit can be obtained using a mass function 
which rises as a power- law with x = 1 down to ~ 0.25-M© 
(M§i4 = 8.5) and then drops all the way to the MS detec- 
tion limit with an exponent x = —0.3. Unfortunately the 
lack of data at the bright end of the WFPC2 LF, due to 
saturation, does not allow us to constrain the MF slope 
in the mass interval which is most sensitive to dynam- 
ical effects, although the comparison between the inner 
and outer mass functions at lower masses already suggests 
strong mass segregation. 

4. The dynamical state of M 4 

Addressing the question as to whether the difference be- 
tween the inner and outer MFs described above can be 
quantitatively accounted for by the two-body relaxation 
mechanism requires a simulation of the dynamical struc- 
ture of the cluster. We have used multi-mass King-Michie 
models, constructed with an approach nearly identical to 
that of Gunn & Griffin (4979), as described in Meylan 
(4987, 1988). Each model is characterized by an IMF in 
the form of an exponential as in Eq. (3), with an expo- 
nent x which is not necessarily constant (i.e. x is allowed 
to vary with mass; see details below), and by three struc- 
tural parameters describing respectively the scale radius 
(r c ), the scale velocity (v s ) and the central value of the di- 
mcnsionless gravitational potential W . We have assumed 
complete isotropy in the velocity distribution. 

Because King-Michie modeling provides a "snapshot" 
of the current dynamical status of the cluster rather than 
its evolution, it is useful to define the mass distribution of 



cluster stars at present (global MF; GMF) as the MF that 
the cluster would have simply as a result of stellar evolu- 
tion (i.e. ignoring any local modifications induced by in- 
ternal dynamics and/or the interaction with the Galactic 
tidal field). Clearly, the IMF and GMF of MS (unevolved) 
stars is the same. For practical purposes, the GMF has 
been divided into sixteen different mass classes, covering 
MS stars, white dwarfs, and heavy remnants. All stars 
lighter than 0.8M Q have been considered still on their 
MS, while heavier stars with initial masses in the range 
[8-100] M have been assigned a final (i.e. post-evolution) 
mass of 1.4 M . White dwarfs have been divided into three 
classes, according to their initial masses, following the pre- 
scriptions of Weidemann (1987) and Bragaglia, Renzini, &: 
Bergeron (1995). Although the analytical form we adopted 
for the IMF would allow the exponent x of the GMF to 
vary, we have initially restricted our investigation to the 
single exponent case. 

From the parameter space defined in this way, we have 
extracted only those models that simultaneously fit both 
the observed surface brightness (SBP) and velocity disper- 
sion (VDP) profiles of the cluster. The fit to the SBP and 
VDP, however, can only constrain r c , v s , and W Ql while 
still allowing the IMF to take on a variety of shapes. It 
is precisely to attempt to break this degeneracy that we 
impose the condition that the model GMF agree with the 
observed LFs. We note here that the local MF at any 
given place inside the cluster originates from the GMF 
and stems directly from imposing the condition that stars 
of different masses be in thermal equilibrium with each 
other. 

In our simulations, we have used the VDP determined 
by Peterson, Rees, & Cudworth (1995). Concerning the 
SBP, we have compared the profile measured by Kron, 
Hewitt, & Wasserman (1984) with that of Tragcr, King, 
& Djorgorvski (4995) and found some discrepancies. We 
have, therefore, adopted an average SBP, defined as the 
parameterized King-type profile (King 1962) that corre- 
sponds to the canonical values of core radius r c = 50" and 
concentration c = 1.6 (Harris 1996; note, however, that 
values as large as r c ~ 90" and c ~ 2 are found in the lit- 
erature, as discussed below). The canonical SBP defined 
in this way gives a reasonable fit to both measurements, 
and we have used the least squares between the canoni- 
cal profile and the observations of Kron et al. (1984) and 
Trager, King, & Djorgorvski (1995) as an estimate of the 
errors. 

Although several models are able to reproduce the ob- 
served radial profiles (SBP and VDP) reasonably well, no 
choice of the exponent x gives local MFs that, converted 
into LFs using the M-L relation described above, agree 
with the data shown in Fig. 14. In fact, our simulations 
show that an IMF like the one given in Eq. (3) cannot 
simultaneously fit the LFs at 2' and 6' radial distance, 
regardless of the value of the exponent x, although each 
LF can individually be approximated with a power-law 
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Fig. 15. The WFPC2 2'(triangles) and 6'(squares) LFs, 
the NICMOS LF (diamonds) in the F814W band, com- 
pared with the theoretical LFs at the same distances, ob- 
tained by the mass stratification of the dynamical model 
described in the text. The slope of the main sequence 
mass segment of the GMF is x = 0.4. The core ra- 
dius and the concentration of the parameterized SBP are 
frmc = 50" and C— 1.6 



IMF with index x ~ —0.8 for the LFs at 2' radius and 
x ~ 0.4 at 6'. (In the latter case, however, the fit is not 
good at magnitudes fainter than M 8 i 4 ~ 9). As an ex- 
ample, in Fig. 15 we show the theoretical LFs (solid lines) 
corresponding to the model that best fits the observations 
at 6': at bright magnitudes (i.e. for masses > 0.3 M Q ) the 
model LF systematically deviates from the observational 
data at 2' by far more than the experimental errors. 

A single-exponent power-law distribution, however, is 
only a rough approximation to the IMF, particularly at 
the low-mass end (see Scalo 1998). And indeed, the deep- 
est LFs today available for NGC 6397 (Paresce, De Marchi, 
& Romaniello 1995; King et al. 1998) and NGC 6656 (Dc 
Marchi & Paresce 1997) suggest that the MF could flatten 
out below ~ 0.25 M Q and possibly drop (in the logarith- 
mic plane). We have, therefore, run some simulations with 
an IMF that follows an exponential law with two values 
of the exponent: x is positive for masses m > 0.25 M Q and 
zero below that mass. Not even in this case, however, does 
the predicted shape of the MF at ~ 2' and ~ 6' from the 
center match the observations simultaneously, regardless 
of the value of x for masses m < 0.25 M Q . 

Although we could explore other possible shapes of 
the IMF, it is unlikely that its form can largely deviate 
from some sort of a slowly varying exponential. Thus, if 
we judge that King-Michie modeling is a viable way to 
describe the dynamical state of M 4, we have to conclude 
that some of the constraints that we are forcing into the 
models are not correct, i.e. that the SBP or the VDP (or 



Fig. 16. The same as Fig. 15 when adopting as GMF 
a power-law with exponents x = 0.2, in the range 
0.8-0.25 M Q , and x = -0.4 for masses < 0.25 M Q 



both) are not compatible with the LFs that we measured 
and which seem rather robust (see Sect. 3). 

We have, therefore, investigated the origin of the dis- 
crepancy between the observed LFs and the local MFs 
produced by the models starting from plots similar to that 
shown in Fig. 15. We have noticed that, if we begin by fix- 
ing the exponent (exponents) of the IMF in such a way 
that the local MF agrees with our measurements at 6' from 
the center, we are tempted to conclude that the scale ra- 
dius of our models is too small. Fixing the exponent in 
that way is a reasonable assumption, since both theoret- 
ical and semiempirical arguments suggest that near the 
cluster half-mass radius (~ half-light ~ 3.7'; Harris 1996), 
the local MF should only marginally differ from the IMF 
(Richer et al. 1991; Vesperini & Heggie 1997). 

In fact, in order for the model LF to reproduce the 
LFs measured at 2', the former has to decrease monoton- 
ically with luminosity, and this is precisely what has been 
shown to take place in the innermost cluster regions (see 
Paresce, De Marchi, & Jedrzejewski 1995; King, Sosin, & 
Cool 1995; De Marchi & Paresce 1996). In other words, 
we would need a scale radius (r c ) sufficiently large so as to 
place the region at r ~ 2' still close enough to the center 
that the LF of its stellar population is compatible with an 
inverted function. The canonical value of r c = 50" is not 
large enough for this to happen. 

Through repeated trials, we have found that by forc- 
ing the core radius r c to take on the value of ~ 90" (a fig- 
ure still compatible with the measurements of Peterson & 
King (1975) suggesting r c ~ 93" from the fit to their SBP) 
both sets of LFs can be simultaneously fit by an exponen- 
tial IMF with index x ~ 0.2 in the range 0.8-0.25 M Q and 
which drops with x ~ —0.4 at lower masses (Fig. 16). The 
values of the concentration ratio (c ~ 1.5), the total cluster 
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Fig. 17. The same as Fig. 15 when the radial LFs are fit- 
ting a SBP with r c = 92" and c = 1.9. The adopted GMF 
has main sequence slopes x = 1, in the range 0.8-0.25 M Q , 
and x = —0.4 for masses < 0.25 M Q 



mass(~ 8.8 x 1O 4 M ), and the index of the IMF above 
0.8 Mq (x ~ 1.4) are also reasonable. In the mass range 
0.2 — 0.6 M , the slope of the IMF suggested by this model 
is shallower than the x ~ 1, value proposed by De Marchi 
& Parcscc (1997) for NGC 6397 and NGC 6656. Neverthe- 
less, it would still be possible to find a model (Fig. 17) that 
uses the IMF of De Marchi & Paresce (1997) and simul- 
taneously fits the LFs observed at ~ 2' and ~ 6', and this 
would require the core radius to grow to <~ 90" and the 
concentration parameter to take on the value of c ~ 1.9 
(both values still compatible with the data of Peterson 
& King (1975) and Kron ct al. (1984), respectively). The 
difference between the IMFs used in Figs. 16 and 17 and 
between the set of parameters defining our models can be 
taken as an indication of the uncertainties affecting our 
approach, and can be traced back to the indetermination 
accompanying the SBP and the ensuing value of r c . It also 
shows how difficult it is to infer the shape of the GMF 
starting from the observations when the knowledge of the 
structural parameters of the cluster is so uncertain. 

We would also like to point out that the results given 
above and shown in Figs. 16 and 17 should anyhow be 
taken with care: they stem from the assumption that King- 
Michie modeling with the limitations that we applied to 
it is adequate to describe the dynamical state of M 4, and 
therefore the good fits that we obtain do not necessarily 
guarantee that these models reflect the physical state of 
the cluster. And indeed, from the kinematics of its orbit 
(Dauphole et al. 1996) one would argue that M4 should 
have undergone extensive tidal stripping and shocks since 
it penetrates very deeply into the Galactic bulge and its 
orbital plane is one of the closest to the Galactic disk. 
These effects might have strongly modified the local mass 



distribution (particularly in the outer regions) and the 
two-body relaxation mechanism that we adopted could 
no longer apply. We are grateful to Pavel Kroupa, the 
referee, for having pointed out that, among the possible 
problems with the assumptions underlying our modeling, 
there is the assumed complete isotropy in the velocity dis- 
tribution of the cluster members. Stellar-dynamical mod- 
cling of isolated globular clusters shows that the velocity 
dispersion becomes anisotropic as a function of time (see 
Fig. 4e in Spurzem & Aarseth 1996). Models that include 
a tidal field, however, suggest that the anisotropy is re- 
duced because of rapid loss of stars having radial orbits 
(Takahashi, Lee, & Inagaki 1997). Nevertheless, it is not 
clear if the tidal field is realistically taken account of in the 
Takahashi et al. (1997) Fokker-Planck models, since in the 
above mentioned Fig. 4e of Spurzem & Aarseth (1996) it 
can be seen that the anisotropy would also be pronounced 
well inside the tidal radius. Moreover, GCs may be born 
with a significant anisotropic velocity dispersion (Aarseth, 
Lin, & Papaloizu 1988). 

With this caveat in mind, we can quantitatively con- 
clude that the LFs observed at ~ 2' and ~ 6' away from 
the center of M 4 can be traced back to the same GMF, 
which appears to have been locally modified as a result 
of mass segregation due to two body relaxation. The un- 
derlying IMF is most likely represented by a slowly rising 
exponential (x ~ 0.2-1) that flattens out and possibly 
drops below ~ 0.25 M Q . Any attempt to better define the 
shape of the IMF and to characterize more thoroughly 
the relaxation mechanism would require a more precise 
measurement of the surface brightness (or, better, radial 
density) profile and the observation of the LF of MS stars 
at four or five locations inside the cluster, so as to simul- 
taneously constrain both the core radius (for a given mass 
class) and the mass spectrum of all luminous objects as a 
function of the distance from the center. This project can 
be efficiently and effectively carried out with a wide field 
imager at a 10-m class telescope like the VLT. 

5. Summary 

Our main conclusions can be summarized as follows: 

1 . We have measured the visible and infrared LFs of M 4 
at 2'and 6' from the cluster center down to I ~ 10 with 
photometric completeness always better than 50%. 

2. The LFs measured at 2' with the two instruments agree 
very well with each other and show a slow decrease 
with luminosity. 

3. The LF measured with the WFPC2 at 6' radial dis- 
tance shows an increase with decreasing luminosity up 
to a peak at I ~ 9 followed by a drop all the way to 
the detection limit. 

4. The difference between the two sets of LFs is qual- 
itatively compatible with the effects of mass segre- 
gation. It can, however, be quantitatively reproduced 
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using a King-Michie model only when adopting clus- 
ter structural parameters which, although compatible 
with those available in the literature, deviate from 
their canonical value. The IMF required for the model 
to fit has an exponential slope of ~ 0.2-1.0 in logarith- 
mic mass units between 0.2 and 0.8 M Q . 

5. Our observations require that the IMF determined in 
this way flatten out and possibly drop below 0.2 M . 

6. For the first time an infrared WD sequence has been 
unambiguously identified in a globular cluster. The 
WD locus appears to be a color sequence instead of 
a luminosity sequence as in the optical. 
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